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Abstract 

We study the effective string corrections to the inter-quark potential at finite 
temperature by simulating the SU (2) lattice gauge theory in four dimensions. 
We provide the first numerical evidence that the logarithmic correction to 
the potential, which was recently proposed to be a signature of the effective 
string at finite temperature, is universal also in (3 + 1)D gauge theory, thus 
extending previous results limited to the (2 + 1)D case. 
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1. Introduction 

A precise understanding of the mechanism responsible for color confine- 
ment is still lacking, however there is by now little doubt that the interaction 
between a static quark- ant iquark pair can be described by an effective string 
theory. From this point of view, initiated by the seminal studies [H-El], the 
presence of a flux tube joining the static quark-antiquark pair is given for 
granted and its physical consequences are examined. 

At zero temperature the best known consequence of the effective string 
description is the presence of the so called "Luscher term" [l|-[3| , i. e. the pres- 
ence in the long distance inter-quark potential of a term proportional to 1/r 
(r being the inter-quark distance), with a known coefficient which depends 
only on the dimension of the space-time (in particular it does not depend 
on the gauge group) and it is not influenced by higher order corrections of 
the effective string action. As a consequence of these general properties, the 
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presence of the Liischer term in the inter-quark potential can be considered 
as a clear signature of the effective string behaviour of the flux tube. Be- 
cause of its theoretical interest much studies have been devoted, in the field 
of lattice gauge theory simulation, to detect the presence of the Liischer term 



The effective string picture can be used to describe the inter-quark po- 
tential also at nonzero temperature T (that must clearly be lower than the 
deconfinement temperature T c ): in this regime the dominant string correc- 
tion to the potential is linear in r and decreases the value of the string 
tension. This is temperature dependent and ultimately vanishes at decon- 
finement. In Q it was recently argued that an universal signature of the 
effective string at finite temperature, analogous to the Liischer term at zero 
temperature, is given by the logarithmic correction to the inter-quark po- 
tential. This conclusion was supported by numerical evidence obtained by 
performing simulations in D = 2 + 1 lattice gauge theories for various gauge 
groups. It was also noted that these numerical simulations are particularly 
time consuming, since the logarithmic correction is washed out by the pro- 
jection to zero transverse momentum of the Polyakov loop, a common trick 
used in order to reduce the statistical errors (see e.g. 0). 

The aim of this paper is to show that the universal form of the logarithmic 
correction to the inter-quark potential at finite temperature (as predicted 
by the effective string picture) is consistent with the numerically generated 
lattice gauge data also for the D = 3 + 1 SU{2) lattice gauge theory. 

2. The lattice setup and the effective string prediction 

We studied the SU(2) Yang-Mills theory discretized on an isotropic cubic 
lattice of spacing a and volume V = a A N t x N% with the usual Wilson action 



Here P^ v {x) is the parallel transport along the elementary loop of the lattice 
(plaquette) oriented along the /iz/-plane and positioned at the lattice site x. 
The sum is extended over all the plaquettes. 

In order to generate an ensemble of configurations with the statistical 
distribution e~ s , with S given by Eq. (p]), the usual Markov chain Monte 
Carlo approach has been used. The updates were performed by using a 
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mixture of heatbath [10 , and overrelaxation [l2[ steps in the ratio of 1 to 
5. 

The Polyakov loop P is defined as the trace of the parallel transport along 
the compactified temporal direction, normalized by the number of colors. The 
inter-quark potential is extracted from the correlators of the Polyakov loop 
according to the relation 

7(r) = ~log<P*(0)P(r)> (2) 

and the correlators of the Polyakov loops are computed by using the multi- 
level algorithm developed in 13| together with the multihit method 



When using the multilevel algorithm, usually the Polyakov loop correla- 
tors at all the needed distances are computed by means of the same ensemble 
of configurations. This introduces strong correlations between the estimated 
values of the potential V(r) calculated at different r's. Because of the limited 
RAM available in the hardware used to perform the simulations, we could 
not adopt this strategy: different runs were used to estimate the Polyakov 
loop correlators at different distances. As a consequence the error analysis 
is simplified, since cross-correlations are not present. Statistical errors are 
evaluated by the bootstrap method (see e.g. [15|). 

In finite temperature gauge theories the form of the inter-quark potential 
predicted by the effective string picture is given, up to an additive normal- 
ization, by the expression (see p, 0]) 

aV{r) =aar + ^-^ log (-) + — + 0(a 3 r" 3 ) (3) 
2N t \aJ r 

Here a is the finite temperature string tension and the second term is the 
logarithmic correction we are interested in. As noted in the introduction the 
coefficient of the logarithmic term is independent of the microscopic details 
of the gauge theory and it is related only to the space-time dimensionality 
D and to the temperature T = l/(aN t ). The term c^a/r is the first non- 
universal correction. 

The expression in Eq. ([3]) has been verified with high accuracy in |7| for 
different D = 2 + 1 lattice gauge theories. Since the D — 2 proportionality 
of the logarithmic term in Eq. is a non-trivial prediction of the effective 
string model, we think it is important to numerically check Eq. also in 
four dimensional lattice gauge theories. 
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3. The numerical results 



In order to test the effective string prediction we used two ensembles 
of data generated from lattices of various temporal extents and different 
coupling values. 

The first ensemble, which will be denoted by / in the following, is obtained 
by using a lattice of size 4 x 40 3 and coupling (3j = 2.1768, which is the 
deconfinement coupling on lattices with temporal extent N t = 3 (see e.g. 



|16I |) and thus corresponds to a temperature Tj = |T C . 

The second ensemble, denoted by II, is obtained by using a 6 x 40 3 lattice 
at fin — 2.2986. This value of the coupling is the critical one on lattices with 
temporal extent N t 



4 (see e.g. [16 



2rp 

3 ±r, 



thus Tn 

The collected statistics was of 10 3 statistically independent measurement 
for each inter-quark distance and for each of the two ensembles. Numerical 
simulations have been performed using GRID resources provided by INFN. 



r/a 


aVi{r) 


aV H {r) 


4 


1.76447(21) 


1.25000(11) 


5 


1.99665(30) 


1.38839(15) 


6 


2.21675(38) 


1.51967(19) 


7 


2.43000(49) 


1.64566(24) 


8 


2.63631(64) 


1.76808(30) 


9 


2.83943(61) 


1.88721(54) 


10 


3.0371(10) 


2.00563(83) 


11 


3.2404(20) 


2.1211(12) 


12 


3.4327(40) 


2.2357(22) 


13 


3.6287(85) 


2.3483(42) 


14 


3.828(19) 


2.4689(78) 



Table 1: Extracted values of the potential V(r) (in lattice units) for the two ensembles. 



The computed values of the inter-quark potential are reported in Tab. (JT]) 
and, for the reason explained in the previous section, all data at different 
values of r/a are completely independent. Because of that, in order to test 
Eq. ([3]), it is more convenient to fit directly the data for the potential instead 
of using its first derivative Q(r) = as is usually done to reduce the amount 
of the cross-correlations. 
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We fitted the values of the potential according to 



r /r\ c 3 a 

aV(r) = c + ci- + c 2 log - H 

a \aJ r 



(4) 



and obtained the best fit parameters reported, together with the reduced x 2 
of the fit, in Tab. 

Since the corrections to the potential predicted by the effective string 
model are given by a power series in l/(cxr 2 ), where a is the finite temperature 
string tension, in order to be consistent with the theory the fitting procedure 
has to be restricted to some region r > r min . In previous studies (starting 
with 17] and up to the recent 0]) it was shown that r min can be safely 
defined by crr^ in = 1.5. We verified that for both the considered ensembles 
the inequality r min < 4a is satisfied; the good agreement of the data with the 
fit function for r > r min is a self consistency check for the chosen r min value. 

We also verified that consistent results are achieved by fitting the deriva- 
tive of the potential instead of the potential itself, however the estimates for 
the parameters obtained in this way are less precise than the ones reported 
in Tab. ©. 

The data for the potential are depicted in Fig. (JT|) together with their best 
fit, while in Fig. ([2]) the derivative of the potential is shown for the ensemble 
II. 





/ 


II 


Co 


0.840(78) 


0.631(41) 


Cl 


0.1763(35) 


0.0997(18) 


c 2 


0.200(43) 


0.165(23) 


c 3 


-0.23(13) 


-0.037(67) 




1.4 


0.42 



Table 2: Best fit parameters for the inter-quark potential obtained by using the functional 
form Eq. (gj). 



The theoretical value predicted for c 2 by the effective string model is 
c l 2 = -f^- For the ensemble / the theoretical expectation is thus cf 1 = 0.25, 
which is consistent with the measured c 2 = 0.200(43). For the ensemble // 
the effective string value is c 2 ft = 0.166 which is in good agreement with the 
lattice estimate c 2 = 0.165(23). 

The parameter ci of the fit is the finite temperature string tension. In 
the Nambu-Goto effective string model this can be explicitly written as a 
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Figure 1: Plot of the potential V(r) (in lattice units) together with the fit for the two 
ensembles. 



function of the zero temperature string tension and it is given by (see 181 ]) 



a 2 a 



a 2 aQi 



2tt 



3N 2 a 2 a 



(5) 



where uq is the zero temperature string tension. By using the values reported 



m 



16] for the string tension at zero temperature for the couplings j3 
and (3 = 2.2986 we get from Eq. (j5J) the estimates 



2.1768 



af G) = 0.1886(84) af^ 



0.1013(14) 



(6) 



which are consistent with the best fit values reported in Tab. (121 . 

As a last observation we notice that the coupling used in the ensemble /, 
Pi = 2.1768, is below the value /3 = 2.2 which is usually reported for the so 
called "bulk" or "roughening" transition. Although for some lattice actions 
this is a real phase transition (see e.g. [19]) for the SU(2) lattice gauge 
theory with the Wilson action it is just a smooth analytical crossover, easily 
seen e.g. in the plaquette susceptibility. This crossover signals the transition 
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Figure 2: Plot of Q(r), the derivative of the potential, (in lattice units) together with the 
fit for the ensemble II. 



from the strongly coupled regime, where lattice artefacts are dominant and 
the strong coupling expansion is convergent, to the weakly coupled phase, 



where the continuum limit is smoothly reached (see e.g. 20] and 10] for the 
first numerical results). 

Our results concerning the ensemble / show that, although it is defined 
only in the continuum, the effective string model is consistent with the be- 
haviour of lattice gauge theories also for values of the coupling where strong 
discretization corrections are to be expected. The fact that the effective string 
description is compatible with the data also below the roughening transition 
is not to be seen as an internal inconsistency. Since the roughening transi- 
tion is not, in the studied theory, a real phase transition, its location is not 
a well defined parameter and the departure of the effective string prediction 
from the lattice data smoothly take place in a quite large interval of coupling 
values. 
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4. Conclusions 

The effective string model is a simple theoretical picture which encloses 
much of the physics of non abelian gauge theories and give testable prediction 
on the behaviour of the inter-quark static potential. 

A clear signature of the effective string at zero temperature is given by 
the well known Liischer term. It is a recent observation that a similar role is 
played, at finite temperature, by the logarithmic correction to the potential. 
This was up to now confirmed by lattice results in gauge theory only in the 
2 + 1-dimensional case. 

We extended these results by providing numerical evidence that the same 
is true also for the SU(2) four dimensional lattice gauge theory. We also 
noticed that the observed behaviour of the potential is well described by the 
effective string model also at very large coupling, where large lattice artefacts 
are expected to be present. 

It can be noted that the data of the ensemble II are in better agreement 
with the theoretical expectations than the one of the ensemble I. This can 
in principle be related to the different temperature (higher for the ensemble 
I) and to the lattice artefacts (larger for the ensemble I). It is indeed known 
that the Nambu-Goto string description is not valid near the deconfinement 
transition (see js| for a study of the string tension in (2 + \)D SU{2) near 
the deconfinement) and it is clear that lattice artefacts are not taken into 
account in the effective string picture. Further studies will be needed to gain 
a better understanding of this point and, more generally, of the limits of 
applicability of the effective string model to lattice gauge theory. 
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